NUMERICAL REGULARIZATION FOR SDEs: CONSTRUCTION OF NONNEGATIVE SOLUTIONS

In the numerical solution of stochastic diierential equations (SDEs) such appearances as sudden, large uctuations (explosions), negative paths or unbounded solutions are sometimes observed in contrast to the qualitative behaviour of the exact solution. To overcome this dilemma we construct regular (bounded) numerical solutions through implicit techniques without discretizing the state space. For discussion and classiication, the notation of life time of numerical solutions is introduced. Thereby the task consists in construction of numerical solutions with lengthened life time up to eternal one. During the exposition we outline the role of implicitness for this`process of numerical regularization'. Boundedness (Nonneg-ativity) of some implicit numerical solutions can be proved at least for a class of linearly bounded models. Balanced implicit methods (BIMs) turn out to be very eecient for this purpose. Furthermore, the local property of conditional positivity of numerical solutions is shown constructively (by special BIMs). The suggested approach also gives some motivation to use BIMs for the construction of numerical solutions for SDEs on bounded manifolds with`natural conditions' on their boundaries. Finally we suggest to apply these methods to population dynamics in Biology, innovation diiusion in Marketing and to mean reverting processes in Finance, such as stochastic interest rates. 1 This contribution exhibits a rst trial to numerical analysis of Stochastic Diierential Algebraic Equations (SDAEs) with nonanticipating algebraic constraints. It is an improved version of discussion paper no. 40 refereed and published within Sonderforschungsbereich 373 at Humboldt University in Berlin (1994).